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Abstract 

In the framework of Brans-Dicke (BD) theory, the first part of the present study determines 
the time dependence of BD parameter, energy density and equation of state (EoS) parameter 
of the cosmic fluid in a universe expanding with acceleration, preceded by a phase of 
deceleration. For this purpose, a scale factor has been chosen such that the deceleration 
parameter, obtained from it, shows a signature flip with time. Considering the dark energy to 
be responsible for the entire pressure, the time evolution of energy parameters for matter and 
dark energy and the EoS parameter for dark energy have been determined. An effective 
interaction term, between matter and dark energy, has been proposed and calculated. Its 
negative value at the present time indicates conversion of matter into dark energy. Using this 
term, the time dependence of the rates of change of matter and dark energy has been 
determined. It is found that the nature of dependence of the scalar field upon the scale factor 
plays a very important role in governing the time evolution of the cosmological quantities 
studied here. The present study provides us with a simple way to determine the time 
evolution of dark energy for a homogeneous and isotropic universe of zero spatial curvature, 
without involving any self-interaction potential or cosmological constant in the formulation. 

Keywords: Brans-Dicke Gravity, Cosmological Parameters, Dark Energy, Matter-Dark 

Energy Interaction Term 

1 Introduction 

The accelerated expansion of the universe is one of the most interesting and important 
phenomena in the field of cosmology that have been obtained through astrophysical 
observations in recent times (Riess et al. 1998; Riess et al. 2001; Bennett et al. 2003; 
Komatsu et al. 2011; Percival et al. 2010). An exotic form of energy, with a negative 
pressure, has been found to be responsible for the accelerated expansion of the universe. This 
energy is known as dark energy (DE). In the fields of physics and astronomy, an extensive 
research is now taking place, throughout the world, on the nature and dynamics of DE. A 
number of models have been proposed to explain this accelerated expansion of the universe, 
following a phase of deceleration. In most of the models, DE is represented by cosmological 
constant A (Carroll 2001). 

The present article is based on Brans-Dicke (BD) theory of gravitation. The BD 
theory is characterized by a scalar field ( <p ) and a dimensionless coupling parameter (m) that 
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govern the dynamics of space-time geometry. It can be regarded as a natural extension of the 
general theory of relativity which is obtained in the limit of an infinite oo and a constant value 
of <p (Sahoo and Singh 2002). The BD theory of gravity can be regarded as one of the most 
important theories, among all prevalent alternative theories of gravitation, which have very 
successfully explained the early and late time behaviours of the universe and solved the 
problems of inflation (Linde 1990). As an extension of the original BD theory, a generalized 
version was proposed, where a) is regarded as a function of the scalar field 0 (Bergmann 
1968; Wagoner 1970; Nordtvedt Jr. 1970). Several models regarding the expanding universe 
have been formulated on the basis of this theory (Sahoo and Singh 2002; Banerjee and 
Ganguly 2009; Chakraborty and Debnath 2009; Jamil and Momeni 2011; Satish and 
Venkateswarlu 2014). 

In the present study, BD field equations (for a spatially flat, homogeneous and 
isotropic space-time) have been used to determine time dependence of the coupling parameter 
(o>), equation of state (EoS) parameter (y) and energy density (p). Time dependence of a> was 
previously studied by many groups, using scale factors that have power law dependence upon 
time (Sahoo and Singh 2002; Jamil and Momeni 2011; Satish and Venkateswarlu 2014). 
These scale factors lead to time independent deceleration parameters that are not consistent 
with the interpretations of recent astrophysical observations. According to astrophysical 
observations, the universe has made a transition from a phase of decelerated expansion to a 
phase of accelerated expansion, implying that the deceleration parameter has a dependence 
upon time (Riess et al. 2001; Chand, Mishra and Pradhan 2016). Taking this fact into 
consideration, we have used a scale factor (a) that leads to a time dependent deceleration 
parameter which changes sign from positive to its present negative value. For simplicity, we 
have assumed power law dependence of co upon <p and also of <p upon a. Time dependence 
of p and the EoS parameter (y), in the present study, are consistent with the results obtained 
from some recent studies based on general relativity with anisotropic space-time 
(Pradhan, Amirhashchi and Saha 2011; Yadav, Rahaman and Ray 2011). Two theoretical 
models have been proposed to determine the time variation of density parameters (fi) for 
matter and dark energy. The results from this model are in agreement with those obtained 
from studies based on completely different premises (Das and Al Mamon 2014; Paul, Thakur 
and Ghose 2010). An effective interaction term ( Q e ff) between matter and dark energy has 
been calculated here and its time dependence has been studied. Its time variation is found to 
be consistent with the findings of some other studies (Zadeh, Sheykhi and Moradpour 2017; 
Cueva and Nucamendi 2010). In terms of Q e // we have determined the time rate of change of 
densities for matter and dark energy. 


2 Metric and Field Equations 


The action, in the generalized Brans-Dicke theory, is expressed as, 

a = - + + (i) 


In equation (1), g is the determinant of the metric tensor g^ v , L m is the Lagrangian for matter, 
R denotes the Ricci scalar, cp is the Brans-Dicke scalar field and w is a dimensionless 
coupling parameter. This parameter is regarded as a function of <p in generalized BD theory. 
The variation of action in equation (1) leads to the following field equations. 


R 


fiV 



02 Mv 2 9 nv 4 > ,/?0 ^ 0 0/i;v 



( 2 ) 
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£!|-V 2 c|) = ' 

at 2 ^ 2cj+3 


( 3 ) 


R^ v is the Ricci tensor and T fiv is the energy-momentum tensor. In the above equations, a 
semicolon stands for a covariant derivative and a comma denotes an ordinary derivative with 
respect to . 

The energy-momentum tensor (7^ v ) for the cosmic constituents is given by, 

Tfiv = (P + P)u^u v + p^p (4) 


Here, p denotes the energy density, p is the isotropic pressure, u v is the four-velocity vector 
and T is the trace of Tj. This form of the energy-momentum tensor is based on an assumption 
that the total matter-energy content of the expanding universe is a perfect fluid. In a co¬ 
moving coordinate system, it v = (0,0,0,1), having the characteristic of g^ v u^u v = 1. 

The line element for a homogeneous and isotropic universe, in Friedmann-Robertson-Walker 
(FRW) cosmology, can be expressed as, 


ds 2 = —dt 2 + a 2 (t) 


dr + r 2 dd 2 + r 2 sin 2 dd% 2 

. 1 — k.v^ 


( 5 ) 


Where a(t) denotes the scale factor, t is the cosmic time, k is the spatial curvature 
parameter. Three spherical polar coordinates are denoted by r, 6 and respectively. Here k 
describes the closed, flat and open universes corresponding to k = 1, 0,-1 respectively. 

The following three equations are the field equations of BD theory for a universe filled with a 
perfect fluid, in the geometry of FRW space-time expressed by equation (5). These equations 
have been obtained from the field equations (2) and (3) in which T^ v from equation (4) has 
been incorporated and the components of the metric tensor g fiv have been taken from the line 
element expressed by equation (5). 


^ a 2 +k ^ $ _ ^(0) <P 2 _ £ 

a 2 a <p 2 cp 2 <p 


( 6 ) 


2 a a 2 +k u)(<p) <p 2 2 -^ + ^= _ YJL 

a a 2 2 <p 2 a <p <p <p 


( 7 ) 


£ , 3 _ p (i-3y) _ 4> £ 

tp dtp tp 2w + 3 cp 2o)+3 


( 8 ) 


In these field equations, y(— P /p) is the equation of state (EoS) parameter for the cosmic 
fluid, which has been treated as a function of time in the present study. 


3 Theoretical Model 


Combining equations (6), (7) and (8), for zero spatial curvature (k = 0), one obtains, 


ri) + ( 2 | + 62-J)a.-6(4 + 2)f=0 

V cp a (bj \a z aj cp 


( 9 ) 


In the present study we have assumed the following ansatzes for <p and o>. These parameters 
have been assumed to have power-law relations with a and cp respectively. 
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“ = “•(*)" (11) 

The expression of <p in equation (10) has been taken from some earlier studies in this regard 
(Banerjee and Ganguly 2009; Chand, Mishra and Pradhan 2016).The reason for choosing the 
above empirical relation for co is the fact that this is regarded as a function of the scalar field 
in the generalized Brans-Dicke theory (Wagoner 1970; Nordtvedt Jr. 1970). Combining 
equation (10) with (11), one gets. 

Using the expressions (10) and (11) in equation (9), one obtains, 

mcon + (n + 4 — 2 ^ (1 — q) = 0 (13) 


Here, q (= —aa/a 2 ) is the deceleration parameter. 

Writing m = o> 0 and q — q 0 (i.e., their values at t — t 0 ) in equation (13), we get, 

6(1 -q 0 ) n+4-2q 0 

771 —- 

n 2 a) 0 n 

Combining equation (6) with (10) and taking k = 0 (for flat space), one gets, 
p = <pH 2 ^3 + 3n — 

Here, H (= a/a) is the Hubble parameter. 

Replacing all parameters in equation (15) by their values at t — t 0 , one obtains, 

= n 2 H 0 2 (^^° 2 3ntf 0 2 — Pol 1 0o ) 

Using equation (16) in equation (14) one gets, 

3(1 -q 0 )H 0 2 n+4-2q 0 

771 —■ - 

3H o 2 +3nH o 2 -p o /0o n 


(14) 


(15) 


(16) 


(17) 


Equations (16) and (17) show that both m 0 and m depends upon the parameter n, which 
determines the time dependence of 0, as per equation (10). 

Using equations (10) and (15) in (7), for k — 0 (zero spatial curvature) and writing ^ = 
—qH 2 one obtains the following expression of the equation of state parameter. 

_ 2q-l-0.5a>n 2 -n-n 2 +nq /ion 

^ 3 + 3n-0.5a>n 2 1 ' 


The value of the equation of state parameter at the present time is therefore given by, 


To = 


2<7 0 -1-O.5o> 0 n 2 -n-n 2 +nq 0 
3+3n-0.5a) o n 2 


(19) 
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To use the expression of a) (eqn. 11), the values of o> () and m should be taken from equations 
(16) and (17) respectively. The value of m, required for the expressions of p and y. in 
equations (15) and (18) respectively, should then be taken from equation (11). The value of 
0, required for equation (15), should be taken from equation (10). 

To determine the time dependence of several cosmological parameters (<p,a), p,y), following 
empirical expression of the scale factor has been used in the present model. 

a = a 0 £ , xp[a{(t/t 0 )^ — l)] (20) 

This scale factor has been so chosen that it generates a deceleration parameter (q — —da/a 2 ) 
which changes sign from positive to negative, as a function of time. This signature flip 
indicates a transition of the universe from a phase of decelerated expansion to a phase of 
accelerated expansion, in accordance with several recent studies based on astrophysical 
observations (Riess et al. 1998; Riess et al. 2001).We have taken a 0 = 1 for all calculations. 
Here a, /? should have the same sign to ensure an increase of the scale factor with time. Using 
equation (20), the Hubble parameter (H) and the deceleration parameter (q) are obtained as, 



a 


q = - 


a/3 / t\@ 1 

to VtJ 


da 


1 + 


l~l 

a/3 



( 21 ) 

( 22 ) 


For 0 < /? < 1 and a > 0 one finds that, q —> +oo as t -> 0 and q -> — 1 as t oo, showing 
clearly a signature flip of q with time. 

Taking H — H 0 and q — q 0 , at t — t 0 , the values of the constants a and /? are obtained as, 


l - tf 0 t 0 (l+<7 0 ) 

P — 1 — H 0 t 0 ( 1 + <7 0 ) 


(23) 

(24) 


It is often necessary to find the evolution of a cosmological quantity as a function of the 
redshift parameter (z) where z = ^ — 1. Using equation (20) one gets the following 
expression as a relation between redshift parameter and time. 

z = ^-l = Exp[-a{(t/t o y - 1)] - 1 (25) 


Using equations (20-24), one can determine the time dependence of <p, a>, p and y from (10), 
(12), (15) and (18) respectively. 

The values of different cosmological parameters used in this article are: 

72— 

H 0 = = 2.33 x 10~ 18 sec~ 1 , q 0 = -0.55, p 0 = 9.9 x 10~ 27 Kgm~ 3 , H DO =0.7 

(Po — — — 1.498 x 10 w Kgs 2 m~ 3 , t 0 — 1.4 x 10 10 Years — 4.42 x 10 17 s 

G o 


4 Determination of Density Parameters 

The total pressure (P) of the entire matter-energy content of the universe is contributed by 
dark energy because, the whole matter content (dark matter + baryonic matter) is regarded as 
pressureless dust (Banerjee and Ganguly 2009; Farajollahi and Mohamadi 2010). Thus we 
can write, 
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P = yp = y D p D 


(26) 


Here, y D is the EoS parameter for dark energy and p D is the density of dark energy. 
Using equation (26), the density parameter for dark energy is given by, 

P Yd 

Let us assume the following empirical expression for the EoS parameter of dark energy. 


(27) 


Yd = Y fit) (28) 

Here fit) is a function of time for which we have chosen two empirical forms in the 
following two models. The reason for assuming the above ansatz is the remarkable similarity 
of the nature of time evolution of the EoS parameter for dark energy (Vinutha et al. 2016; 
Pasqua and Chattopadhyay 2016) with that for the total matter-energy content 
(Pradhan, Amirhashchi and Saha 2011; Yadav, Rahaman and Ray 2011) of the universe. Here 
fit) is assumed to be a function that connects y D to y. 


4.1 Model -1 


We propose the following form of f(t) for model-1 of density parameter calculation. 

= A (if 

Here A and p are dimensionless constants. 

Using equations (28) and (29) in (27) we get, 


(29) 


fl D = p -° = JL = a - 1 (-) 11 

P Yd \t 0 J 


(30) 


Taking D. d — fl D0 at the present time, i.e., at t — t 0 , in equation (30), we get, 
i 


A = 


^D0 


(31) 


Using equation (31) in (30) we get, 

n D =n D0 {^y (32) 

Using equation (32), the density parameter for matter is given by, 

= p -f = PJ r- = 1 - ^ = 1 - n ^ 0 HV w 

In deriving equation (33), we have used the fact that the present universe is composed mainly 
of matter and dark energy, considering all other forms of energy to be negligibly small (Pal 
2000; Goswami 2017). 
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Using equations (32) and (33), the expressions for densities of dark energy and matter can be 
respectively written as, 

Pd — P^d — P^do (34) 

Pm = P^m = P [l - &D0 {£) (35) 

Using equations (29) and (31) in (28), the equation-of-state parameter for dark energy (y D ) 
becomes, 


= _y_ _ _z_ rty 

fi-DO Vt0/ 


(36) 


Using equation (25) in (32) and (33), the expressions of fl D and £l m can be written in terms of 
redshift (z) as, 


Q d ~ n »o ( 1 + a ln z+J 


-v/P 


n m = i-n D0 (i + i/n^) 


-n/P 


(37) 

(38) 


The present value of fi. D0 is close to 0.7 according to several astrophysical observations (Pal 
2000; Goswami 2017). In the evolution of density parameters of the universe, there was a 
time in the recent past when fl D — fi m = 0.5 and, the corresponding z value was lying 
somewhere in the range of 0 < z < 1, as obtained from some recent studies (Das and A1 
Mamon 2014; Paul, Thakur and Ghose 2010) in this regard. That phase of the universe might 
be same as (or close to) the one when there was a transition from a decelerated expansion to 
an accelerated expansion of the universe because it is the dark energy that is known to be 
responsible for the accelerated expansion. According to a recent study, the transition of the 
universe from a phase of decelerated expansion to its present state of accelerated expansion 
took place in the past at z = 0.6818 which was around 7.2371 x 10 9 years ago (Goswami 
2017). Choosing z c to denote the value of the redshift parameter at which the universe had 
fl D — n m , we get the following expression of /r from equations (37) and (38). 


P = 


(3 ln(2CL D0 ) 

ln\l+-ln(^—)] 
a \z c +lJ\ 


(39) 


Using equations (10), (20), (32) and (33), Pl D and fi. m can be expressed as functions of the 
scalar field ( <p ) by the following two equations. 

a D = n D0 (i + -ml) (40) 

n m = i-n D „(i + -(„£) (4i) 

4.2 Model - 2 


We propose the following form of /(t) for model-2 of density parameter calculation. 
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fit) = A Exp (a£) 

Where, A and A are dimensionless constants. 

Combining equations (27), (28) and (42) we get, 

a ‘ = m m 7 Bxp (- x 3 

Using the fact that 0. D = O. D0 at t — t 0 in (43) we get, 

^ _ Exp(-A) 

^DO 

Substituting for A in equation (43) from equation (44) we get, 
n D = n D 0 Ex P {A( i-£)} 


(42) 


(43) 


(44) 


(45) 


Here fl D0 is approximately equal to 0.7, according to various studies on dark energy (Pal 
2000; Goswami 2017). 

Using equation (45), the density of dark energy is given by, 


Pd = pPl D = pPl D0 Exp {A (l - £)} 


(46) 


Using equation (45), the density parameter for matter ( dark + baryonic ) can be expressed as, 


/-» _ Pm _ P Pd 

n 

P P 


= 1 - n D = 1 - fL D0 Exp {a (l - £)} 


(47) 


Using equation (47), the density of matter can be expressed as, 
Pm = P&m = P [ 1 - UdqExP (l - f)}] 


(48) 


Using equations (28), (42) and (44), the equation-of-state parameter for the dark energy (y D ) 
can be expressed as, 


Combining equations (25), (45) and (47), the expressions of EL D and fl m can be written in 
terms of redshift (z) as, 


fl D = fl D0 Exp 


[ Jr 

ln(z +) 

r A U 

1 a 1 MJ 


n m = i - n D0 Exp [—A {[i - ^12] - 1 } 


(50) 

(51) 


Taking fl D = Pl m at z = z c , one gets the following expression of A from equations (50) and 
(51). 
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A = 




ln(z c + 
a I 


(52) 


Using equations (10) and (20) in (45) and (47), £l D and fi. m can be expressed as functions of 
the scalar field (<p) by the following two equations. 


n D = n D „a P {-A([5^ + i] 1/ ' , - 1 )J 


n m = i - n 


-^h([^ + i] w -i)] 


(53) 

(54) 


4.3 A Comparison of Two Models for Density Parameter Calculation 

For a positive value of z c , both p and A are negative, as obtained from model-1 (eqn. 39) and 
model-2 (eqn. 52) respectively, indicating an increase of Cl D ( and a corresponding decrease 
of H m ) with time. Qualitatively, these two models show the same nature of time dependence 
for H d and H m . with fl D = H O0 and fl m = fl m0 at z = 0 (present time). The time derivatives 
of them, for two models, are respectively given by equations (55) and (56). 

n D =^n D0 (±y ll ~ 1 (55) 

&d = ^n D0 Exp [a (l - £)} (56) 

Since p and A are negative, both expressions of fl D are positive, being proportional to —— at 

to 

the present time. Near the present time, they have almost the same behaviour (fig. 5). 


5 Calculation of an Effective Interaction Term 


The energy conservation equation is expressed as (Sahoo and Singh 2002), 

p + 3 H{p + P) = 0 (57) 

Taking the pressure P — yp — y D p D (eqn. 26) and p — p m + p D in equation (57) one gets, 

Pm + Pd + 3 H[p m + p D ( 1 + y D )] = 0 (58) 

If it is assumed that the two entities, matter and dark energy the universe, have been 
interacting with each other, causing the generation of one of them at the cost of the other, one 
may define a parameter representing their interaction, on the basis of equation (58). The 
interaction term (Q) can be represented by the following equations (Das and A1 Mamon 2014; 
Zadeh, Sheykhi and Moradpour 2017; Farajollahi and Mohamadi 2010). 


(59) 


Pm 3" 3 H Pm — Q 

p D + 3Hp D (l + y D ) = -Q 


(60) 
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A negative value of Q represents a transfer of energy from the matter field to the field of dark 
energy and a positive value of Q represents conversion of dark energy into matter. It is 
evident from equations (59) and (60) that Q has a dependence upon time and, its values, 
obtained from these two equations, would not be the same. Let us denote the values of Q , 
obtained from these two equations, by Q x and Q 2 respectively. 

Using the relations p — p m + p D and p D — p£l D in (59) we get, 

Qi = (1 - n D )p -pti D + 3Hp(l - Q d ) (61) 

y 

Using the relations p D — pil D and y D — — (eqn. 36) in equation (60) we get. 


Q 2 = -a D p - pil D - 3 Hp(VL D + y) 


(62) 


To estimate the difference (AQ) between Qpind Q 2 , we write the following expression. 

AQ = <h-Q 2 =p + 3Hp(l + y) (63) 

To have an estimate of an effective interaction between matter and dark energy, we define a 
parameter Q e // as the average of Qpind Q 2 . Thus we have, 

Qeff = ^ = (j - n D ) p - pil D + S -Hp( 1 - 2H d - y) (64) 


To judge the plausibility of using Q e f f as a measure of interaction between matter and dark 
energy, we have calculated the following expression of A Q/Q e ff, using equations (63), (64). 

A <3 _ p+3H p(l+y) ,sc\ 

/I \ *3 ( OJ ) 

Qeff \^-^>-D)P-P^D + 2 li P^ 1 ~ 2 ^ L D~y) 


Here we have used model-1 of density parameter calculation, because two models predict 
almost the same nature of time dependence of density parameters near the present time (fig. 
5). To determine the time dependence of the interaction term (Q), we have calculated p and 
based on the expressions of p and Vl D in equations (15) and (32) respectively. 

p = p//[n-2(l + ,)- i ^y (66) 


1 



Equations (10) and (11) have been used, along with the definitions of H and q, to calculate p 
(eqn. 66) from (15). Equation (21) has been used to calculate f 1 D (eqn. 67) from (32). 

Using equations (66) and (67), in (61), (62), (64) and (65) one obtains, 

0l = (1 - 1 WPH [n - 2(1 + «,) -+ ££2»(^p + 3Hp(l -fl D ) (68) 
Q 2 = -SI dP H [n - 2(1 + «,) - - 3HP(SI D + y) (69) 
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Qe„ = ()- n D ) pH [n ■- 2(1 + ,) - + 2£0» (^)-o + i Hp(1 - 2 fi D - 7 ) (70) 

If the value of -^- is found to be sufficiently small, it would be quite reasonable to believe 

that Q e ff, defined as an average of Qptnd Q 2 , will truly represent the interaction between 
matter and dark energy. Table-1 shows the values of (Qeff) _ and (A Q/Q e ff) _ for 

t to t to 

several values of the parameter n. Here, the present value of Q e ff is numerically of the order 
of 1(T 44 and it is negative. The present value of A Q/Q e ff is so small (~10 -15 ) that Q e ff, as 
defined by eqn. 64, is likely to describe the interaction with sufficient accuracy. Its negative 
value indicates decay of matter into dark energy at the present time. As n becomes more 
negative, (Q e ff) _ becomes more negative. It implies that a greater rate of change of the 

scalar field ( <p ) has a connection to a larger interaction between mater and dark energy of the 
universe, causing a larger rate of generation of dark energy at the expense of matter. 


TABLE - 1 

Values of Q e // and |AQ/Q e /y | at the present time for several values of the parameter n which 
determines the time variation of the scalar field (0). 

n 

( Qe ff\=t 0 

\Qeff) t=tn 

-1.60 

-1.153 x 10“ 44 

1.296 x 10“ 14 

-1.62 

-1.606 x 10“ 44 

3.719 x 10“ 15 

-1.64 

-2.046 x 10“ 44 

2.433 x 10“ 15 

-1.66 

-2.473 x 10“ 44 

3.020 x 10~ 15 

-1.68 

-2.886 x 10“ 44 

1.021 x 10“ 15 

-1.70 

-3.285 x 10“ 44 

9.549 x 10- 15 

-1.72 

-3.670 x 10“ 44 

4.883 x 10- 15 

-1.74 

-4.042 x 10“ 44 

4.188 x 10- 15 

-1.76 

-4.400 x 10“ 44 

2.602 x 10“ 15 

-1.78 

-4.744 x 10“ 44 

2.001 x 10- 15 

-1.80 

-5.075 x 10“ 44 

1.014 x 10“ 15 

-1.82 

-5.392 x 10“ 44 

6.463 x 10- 15 

-1.84 

-5.696 x 10“ 44 

2.797 x 10“ 15 

-1.86 

-5.986 x 10“ 44 

1.663 x 10“ 15 


Time dependence of p m and p D can be calculated from equations (59) and (60) by taking 
Q = Qeff ■ Use of these two equations for this purpose, instead of (48) and (46), ensures 
consistency with the energy conservation equation given by (57). Thus we can write, 


Pm Qeff 3H Pm 


(72) 


Pd = -Qeff ~ 3Hp D (l + Yd) 


(73) 
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Combining equation (70) with (72) and (73), one gets the following expressions (eqns. 74 and 
75) for the time derivatives of matter and dark energy ( p m and p D ). 


Pm 


= {b n °) pH 


n — 2(1 + q) — 


(jomn° 


6 + 6 n — con 2 


t Ppttp (Ht 0 \ 
t 0 V a/? / 


1 

1^3 


+ ~Hp( 1 - 2- y) - 3Hp(l - n D ) 


Pd — 


1 

2 _n ° 


)pH 


n — 2(1 + q ) 


wmn“ 


6 + 6 n — wn 2 


ppVL D (Ht 0 \ i-/? 


t 0 Va/? 


■ -Hp(l - 2 H d - +) - 3//pH D (1 + y D ) 


(74) 


(75) 


Using these relations we have studied the time dependence of p m and p D graphically (figs. 
13-16). It is evident here that the nature of their time variation must depend upon the 
parameters n and p, where p is connected to the value of z c , as given by equation (39). A 
negative value of p m and a simultaneous positive value of p D imply a conversion of matter 
into dark energy. 


6 Choice of Values for the Parameter n 


The parameter n controls the time evolution of the scalar field (0). Time variations of a>. p 
and y depend on n, as we know that both m and co 0 are functions of n. The scalar field (4>) is 
the reciprocal of the gravitational constant (G) (Banerjee and Ganguly 2009; Roy, 
Chattopadhyay and Pasqua 2013). Hence, using equations (10) and (20), we can write, 

C = £ = t© ” = j-/xp[-na{C/t o y - 1}] (76) 


According to some studies, the gravitational constant increases with time (Pradhan, Saha 
and Rikhvitsky 2015; Saha, Rikhvitsky and Pradhan 2015; Mukhopadhyay, Chakraborty, 
Ray and Usmani 2016). In the expression of G in equation (76), the scale factor (a) is an 
increasing function of time. For negative values of the parameter n, G would be an increasing 
function of time. Therefore, in the present study we have used only negative values of the 
parameter n. 

Using equations (21) and (76), we get, 


G 

G 


—n- = —nH — — n — 

a t 0 



(77) 


Writing equation (77) for the present time, i.e., t — t 0 , one obtains, 

_ _ to_ ( 
a(3 \g) t 


n 


t = t o 


(78) 


Using equation (78), the parameter n can be estimated from the experimental observations 
regarding [6/6] , as reported in articles on time varying gravitational constant (Ray, 

Mukhopadhyay and Dutta Choudhury 2007). According to a study by S. Weinberg, the 


largest possible value of 


t=t 0 


is 4 x 10 10 per year (Weinberg 1972). Using equation (78), 


one gets n -5.5 for this upper limit. 
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TABLE - 2 

Values of co 0 , m , y 0 , y D0 , (-) for several values of the parameter 

yGJ t=t 0 

n which determines the time variation of the scalar field (0). 

n 

m 0 

m 

To 

Ydo 

(z c — 0.7) 

(0 

V y t=t 0 

(Tr _1 ) 

- 1.63 

- 1.514 

- 0.183 

- 1.796 

- 2.566 

1 . 198 X 10“ 10 

- 1.64 

- 1.518 

- 0.168 

- 1.691 

- 2.416 

1 . 205 X 10“ 10 

- 1.65 

- 1.522 

- 0.154 

- 1.588 

- 2.268 

1 . 212 X 10“ 10 

- 1.66 

- 1.525 

- 0.140 

- 1.486 

- 2.123 

1 . 220 X 10“ 10 

- 1.67 

- 1.529 

- 0.127 

- 1.386 

- 1.979 

1 . 227 X 10“ 10 

- 1.68 

- 1.532 

- 0.115 

- 1.287 

- 1.839 

1 . 234 X 10“ 10 

- 1.69 

- 1.535 

- 0.104 

- 1.190 

- 1.700 

1 . 242 X 10“ 10 

- 1.70 

- 1.538 

- 0.093 

- 1.095 

- 1.564 

1 . 249 X 10“ 10 

- 1.71 

- 1.540 

- 0.083 

- 1.001 

- 1.430 

1 . 256 X 10“ 10 

- 1.72 

- 1.543 

- 0.073 

- 0.909 

- 1.299 

1 . 264 X 10“ 10 

- 1.73 

- 1.545 

- 0.064 

- 0.819 

- 1.170 

1.271 x 10“ 10 

- 1.74 

- 1.547 

- 0.055 

- 0.730 

- 1.043 

1 . 279 X 10“ 10 

- 1.75 

- 1.549 

- 0.046 

- 0.643 

- 0.918 

1 . 286 X 10“ 10 

- 1.76 

- 1.551 

- 0.038 

- 0.558 

- 0.796 

1 . 293 X 10“ 10 

- 1.77 

- 1.553 

- 0.031 

- 0.474 

- 0.677 

1 . 301 X 10“ 10 

- 1.78 

- 1.554 

- 0.024 

- 0.392 

- 0.559 

1 . 308 X 10“ 10 

- 1.79 

- 1.555 

- 0.017 

- 0.311 

- 0.444 

1 . 315 X 10“ 10 

- 1.80 

- 1.557 

- 0.011 

- 0.232 

- 0.332 

1 . 323 X 10“ 10 

- 1.81 

- 1.558 

- 0.005 

- 0.155 

- 0.221 

1 . 330 X 10“ 10 

- 1.82 

- 1.559 

0.001 

- 0.079 

- 0.113 

1 . 337 X 10“ 10 

- 1.83 

- 1.560 

0.006 

- 0.005 

- 0.008 

1 . 345 X 10“ 10 

- 1.84 

- 1.561 

0.012 

0.067 

0.095 

1 . 352 X 10“ 10 

- 1.85 

- 1.561 

0.016 

0.137 

0.196 

1 . 359 X 10“ 10 


Table-2 shows the values of 6J 0 , m, y 0 , y D0 , (-) for different values of n which governs the 

w t=t 0 

time variation of the scalar field (0), as per equation (10). For the values of n over the range of 
—1.76 < n < —1.67, the values of the present EoS parameter (y 0 ) in this table are consistent 
with the range of values (approximately —1.7 < y 0 < —0.6) obtained from astrophysical 
observations, as reported in some recent articles on the time variation of EoS parameter, based on 
an anisotropic space-time (Pradhan, Amirhashchi and Saha 2011; Yadav, Rahaman and Ray 
2011). Over this range, the values of m are negative, indicating that co would be less negative as 
0 increases, as evident from equation (11), where co 0 is negative. For n < —1.81 in this table we 
have m > 0, implying that a> would be more negative as 0 increases. Both of these natures of 
variations of m(0) are reflected in the plots of some other studies on Brans-Dicke theory 
(Chakraborty and Debnath 2009; Sharif and Waheed 2012). But the second type of variation is 
not physically valid, since it corresponds to those values of EoS parameter which are not in 

conformity with the astrophysical observations in this regard. The values of (-) in this table 

VG/ t=t 0 

are absolutely within the range of permissible values shown by S. Weinberg (Weinberg 1972). 
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parameter versus time for three different 
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Fig. 10 Plots of effective interaction term 
Q e ff versus time for three different values 
of the parameter n 
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Fig. 11 Plots of effective interaction term 
Q e ff versus redshift (z) for three different 
values of z c 
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Fig. 12 Plots of effective interaction term 
Qeff versus redshift (z) for three different 
values of the parameter n 
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Fig. 13 Plots of p m versus time for three 
different values of z c 
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Fig. 14 Plots of p D versus time for three 
different values of z c 
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Fig. 15 Plots of p m versus time for three 
different values of the parameter n 


:—>—i—>—i—■—r 

—•—i— 1 —i—■—: 

-n = -1.8 . 

\ 

- n = -1.7 


- - - n = -1.6 

* 

z =0.1 ■ 

- 


V 




_1_1_1_1_1_L 

- 


0.7 0.8 0.9 1.0 1.1 1.2 1.3 

t/t„ 


Fig. 16 Plots of p D versus time for three 
different values of the parameter n 
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7 Results and Discussion 

Figure 1 shows the variation of the Brans-Dicke parameter (n>) as a function of time, for 
different values of the parameter n. Here, co has a negative value and it becomes more 
negative with time. For more negative values of n, it decreases less rapidly. This nature of co, 
which is gradually becoming more negative with time, is also obtained from other studies on 
Brans-Dicke theory (Chakraborty and Debnath 2009; Sharif and Waheed 2012). The time 
dependence of co, obtained from the present study, is based on a scale factor which leads to a 
time dependent deceleration parameter (that shows a signature flip with time from positive to 
negative), unlike some recent studies (Sahoo and Singh 2002; Jamil and Momeni 2011; 
Satish and Venkateswarlu 2014) that are based on time independent deceleration parameters. 
The value of co 0 (close to —1.55) and its gradual decrease with time are consistent with the 
results provided in the section of ‘conclusions’ of a recently published article based on 
Kantowski-Sachs space-time (Satish and Venkateswarlu 2014). 

Figure 2 shows the time variation of energy density (p) with time, for different values of the 
parameter n. Here, p rises steeply to a peak and then decreases at a slower rate. More 
negative values of n cause quicker attainment of the peak and also a larger peak value of 
energy density. It falls more steeply for more negative values of n. This behaviour is 
consistent with studies based on General Relativity (Pradhan, Amirhashchi and Saha 2011). 

Figure 3 shows the time dependence of the equation of state (EoS) parameter (y) for three 
different values of the parameter n. These curves show a steep rise initially and then they 
gradually become asymptotic to a small negative value. Saturation is reached faster for less 
negative values of n. The value of y at the present time (t = t 0 ) is —0.73. The ranges of y 
from observational results based on SN la data and SN la data collaborated with galaxy 
clustering statistics and CMBR anisotropy are -1.67 < y < -0.62 and -1.33 < y < -0.79 
respectively. (Knop et al. 2003; Tegmark et al. 2004). These curves are similar to those 
obtained from dark energy models based on LRS Bianchi type-V metric in the framework of 
Einstein’s general theory of relativity (Yadav, Rahaman and Ray 2011). 


Figure 4 shows a single plot EoS parameter, for n — —1.74, on a longer time scale. It shows 
a pattern that is very much similar to the one depicted in one of the earlier studies based on 
dark energy models in anisotropic Bianchi type-I (B-I) space-time, in the framework of 
Einstein’s general theory of relativity (Pradhan, Amirhashchi and Saha 2011). 

The plots in the figures 2-4 show very clearly that the time evolution of energy density and 
EoS parameter, obtained here for FRW space-time in the framework of Brans-Dicke theory, 
are quite in agreement with the results obtained from anisotropic space-time in the framework 
of general relativity. 

Figure 5 shows the variations of density parameters, for dark energy and matter, as functions 
of the redshift parameter, using models-1, 2 (Sec. 4.1, 4.2). Here, the value of /r has been 
taken to be —0.557, causing z c = 0.7, implying that the dark energy took over at around 
z = 0.7. These plots are quite similar in nature to the plots obtained from some recent studies 
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(Das and A1 Mamon 2014; Paul, Thakur and Ghose 2010) which were carried out on 
premises completely different from the present study. 

Figure 6 shows the plots of density parameters, for dark energy and matter, as functions of 
time (Sec. 4.1), for three different values of z c . For higher values of z c (implying less 
negative values of /r, as per eqn. 39), |fl| gradually decreases with time. 

Figure 7 shows the time variation of the EoS parameter for dark energy ( y D ) for three 
different values of the parameter n. These plots are very much similar to the ones obtained 
from studies based on anisotropic space-time in the framework of general relativity 
(Pradhan, Amirhashchi and Saha 2011; Yadav, Rahaman and Ray 2011). Plots with less 
negative values of n take longer time to attain saturation. 

Figure 8 shows the time variations of EoS parameters for total energy and dark energy. Here 
y D has a more negative value compared with y. Their difference decreases with time, as the 
proportion of dark energy increases in the universe. For n = —1.74 the values of y 0 and y D0 
are respectively —0.73 and —1.04 respectively. These values are quite consistent with the 
ranges of values obtained from several astrophysical observations and described in some 
studies regarding dark energy and time varying EoS parameter (Pradhan, Amirhashchi and 
Saha 2011; Yadav, Rahaman and Ray 2011; Das and A1 Mamon 2014). 

The parameter n, in this model, plays an important role in governing the time evolution of 
various cosmological quantities. Both m 0 and m, which controls the value of n>, are 
dependent upon n (eqns. 16, 17). Larger negative value of n means faster change of <p 

(= 1/(7) with time. Taking n — —1.74 we get, = 1.28 x 10 -10 Tr _1 from equation 

\ G -'t=t 0 

(77). This is well within its upper limit, i.e. 4 x 10 -10 Tr -1 , predicted by S. Weinberg 
(Weinberg 1972). 

Figure 9 shows the variation of the effective interaction term Q e ff as a function of time for 
three different values of z c . It is found to decrease with time from an initial positive value, 
gradually becoming negative at a later time. This behaviour of signature flip is consistent 
with the results of a recent study on holographic dark energy (Zadeh, Sheykhi and Moradpour 
2017). For smaller values of z c (which correspond to more negative values of p), Q e f f attains 
its negative phase faster. Its positive phase indicates that dark energy initially decayed into 
matter and, in the later stage, where Q e ff < 0, matter decays into dark energy. 

Figure 10 shows the variation of Q e // as a function of time for three different values of the 
parameter n. Two of these curves decrease with time and change their signs from positive to 
negative. These two curves show that the transfer of energy, between matter and dark energy, 
reverses its direction, producing dark energy at the expense of matter at the later stage. The 
third curve, with n — —1.8, is entirely in the negative domain, becoming gradually closer to 
zero with time. It clearly shows that, there exists a value of the parameter n, below which we 
have an unidirectional energy transfer, from the matter field to the field of dark energy, at a 
gradually decreasing rate. 

Figure 11 depicts the change of Qeff as a function of z, for three different values of z c . It is 
found that, for smaller values of z c , Q e // attains its negative phase faster. 
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Figure 12 shows the variation of Q e ff as a function of the redshift parameter (z) for three 
different values of the parameter n. Like figure 10, one of these three plots also indicate the 
existence of a threshold value for n. below which the direction of energy transfer is always 
from the matter field to the field of dark energy with Q e f j becoming less negative with time. 
Above this threshold of n, the interaction term shows a signature flip from positive to 
negative as time progresses. The curve for n — —1.8 in this figure is similar to the behaviour 
obtained from a study on the interaction term between dark matter and dark energy (Cueva 
and Nucamendi 2010). Since around 80% of the matter content is dark matter, the 
interactions of dark energy with matter consist mainly of the interactions of dark energy with 
dark matter. 

Figure 13 shows the variation of p m as a function of time for three different values of z c . 
Each graph shows a decrease with time, becoming positive to negative prior to the present 
time (t = t 0 ). Smaller values of z c cause faster transition from positive to negative values. In 
an interacting model of matter and dark energy, a negative value of p m means a transfer of 
energy from the sector of matter to that of dark energy. 

Figure 14 depicts the variation of p D as a function of time for three different values of z c . 
Each plot shows a decrease of p D with time. For the larger values of z c , we find a faster fall 
of these graphs. For a conversion matter into dark energy, at the present time, it is necessary 
that p D is positive when p m is negative. Figures 13, 14 show that it is possible to find a value 
of z c for which this condition is satisfied during a span close to the present time (t = t 0 ). 

Figure 15: shows the variation of p m as a function of time for three different values of the 
parameter n. For less negative values of this parameter, these curves show a decrease with 
time, from positive to the negative region. For more negative values of n, p m remains 
negative throughout, with a gradually decreasing slope. For less negative values of n, p m 
becomes negative prior to the present time (t = t 0 ). It means that the matter content started 
decaying into dark energy somewhere in the recent past, causing an accelerated expansion of 
the universe. 

Figure 16 shows the time variation of p D for three different values of the parameter n. For 
z c — 0.1, p D is positive for all n values here (implying a creation of dark energy) and it 
decreases with time. The curves with more negative values of n are lying below the ones with 
less negative values of the same parameter. It is evident that the more negative values of n 
cause slower rate of production of dark energy. 

8 Concluding Remarks 

The present study is based on a spatially flat, homogeneous and isotropic space-time, in the 
framework of Brans-Dicke theory. Here, the time evolution of Brans-Dicke parameter (oo), 
energy density (p) and EoS parameter (y) have been determined from the field equations. The 
choice of the scale factor, that ensures a signature flip of the deceleration parameter with 
time, is likely to give more credibility to the results of this model. Without involving any self 
interaction potential (V4) and cosmological constant (A), the time evolutions of density 
parameters for matter and dark energy have been determined and these results are sufficiently 
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consistent with those obtained through different models and methods (Das and A1 Mamon 
2014; Paul, Thakur and Ghose 2010). The parameters p and A, in model-1 and model-2 (sec. 
4.1, 4.2), govern the time dependence of Q. D , D. m and y D . It has been shown that p and A can 
be determined if one gets information from astrophysical observations regarding the time or 
redshift (z) at which we had 11 D = 11 m in the recent past of the expanding universe. Assuming 
an interaction to be taking place between matter and dark energy, an effective interaction 
term has been calculated on the basis of the energy conservation equation. For some 
combinations of n and z c (or n and p), this term changes its sign from positive to negative, 
showing a change of direction of energy transfer between the sectors of matter and dark 
energy. Its negative value at the present time implies a decay of matter into dark energy 
which is responsible for the accelerated expansion of the universe. The behaviour of this 
interaction term is consistent with the results of some recent studies (Zadeh, Sheykhi and 
Moradpour 2017; Cueva and Nucamendi 2010) carried out in a completely different way. 
Using the expression of Q e jf we have determined the time variations of p m and p D . The 
conditions under which they respectively have negative and positive signs during a span of 
cosmic time, matter may be regarded as decaying into dark energy in that phase of expansion 
of the universe. The parameter n, which actually controls the rate of change of gravitational 
constant with time (as per eqn. 76), determines the nature of time variation of some 
cosmological quantities (co, p and y) mentioned above. It implies that the dependence of the 
scalar field (0 = 1 /G) upon the scale factor (a) (eqn. 10) plays a very important role in 
cosmic expansion. The dependence of co upon (p is controlled by m 0 and m which are all 
functions of n. Positive values of n leads to large negative values of y D(h contrary to its range 
of values close to —1, obtained from recent astrophysical observations (Das and A1 Mamon 
2014). This result is in favour of using negative values of n, implying an increase of 
gravitational constant with time, as per equation (76). This behaviour of gravitational 
constant is also obtained from other studies (Roy, Chattopadhyay and Pasqua 2013; 
Pradhan, Saha and Rikhvitsky 2015). Keeping in view all these facts, one may think of 
improving this model by choosing a relation, between 0 and a, which is different from the 
ansatz of equation (10), in order to have more than one tunable parameter like n controlling 
the time variation of the scalar field (0). As a future project, one may think of determining 
the time dependence of p m and p D from equations (59) and (60) by assuming Q to have a 
form like Q — kH where k is a constant and H is the Hubble parameter. In the second 
equation, y D can be replaced by y D 0 , considering a very slow time variation of y D , as evident 
from figure 8, over a span of time close to t = t 0 . The time dependence of H can be taken 
from equation (21). Using the solutions of these differential equations, (59) and (60), one can 
then calculate II D and II m . taking the energy density (p) from equation (15). It would be 
possible to make an estimate of k by using the known values of fl D0 and Q m0 at t — t 0 in the 
expressions of the density parameters for matter and dark energy. Using equations (10) and 
(19), one would be able to express the density parameters II D and II m in terms of the scalar 
field (0). 
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